Electronic transport in graphene-based structures: an effective cross section approach 
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We show that transport in low-dimensional carbon structures with finite concentrations of scat- 
terers can be modeled by utilising scaling theory and effective cross sections. Our reults are based 
on large scale numerical simulations of carbon nanotubes and graphene nanoribbons, using a tight- 
binding model with parameters obtained from first principles electronic structure calculations. As 
shown by a comprehensive statistical analysis, the scattering cross sections can be used to esti- 
mate the conductance of a quasi-lD system both in the Ohmic and localized regimes. They can 
be computed with good accuracy from the transmission functions of single defects, greatly reducing 
the computational cost and paving the way towards using first principles methods to evaluate the 
conductance of mesoscopic systems, consisting of millions of atoms. 

PACS numbers: 72.10.-d, 72.80.Vp, 73.23.-b 
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Graphene, an effectively two-dimensional material con- 
sisting of a single sheet of carbon atoms, is regarded to be 
a potential candidate for a wide range of future electronic 
devices [H . In order to characterize phenomena affecting 
charge carrier transport in graphene-based systems, ef- 
fective computational methods are required. Particularly 
important is the study of effects that induce a transport 
gap, turning graphene into a semiconductor. 

The mechanisms behind experimentally measured 
transport gaps in graphene nanoribbons (GNRs) are cur- 
rently actively debated. On one hand, it has been sug- 
gested that Coulomb blockade may simificantly limit the 
conductance close to the Dirac point 5| , through bar- 
riers formed by either edge roughness [3| or charged im- 
purities 0, although their role in limiting mobility has 
been questioned [6|. On the other handj^as the phase 
coherence length is very long in graphene [T] , also Ander- 
son localization (AL) may induce a transport gap [S- 13|. 
This is supported by recent experimental results, which 
show a resistance growing exponentially with length in 
some GNRs 
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The low-energy band structure of graphene has two 
nonequivalent valleys, and due to the linear dispersion 
relation, intravalley scatte ring will result in antilocaliza- 
tion instead of localization 15|, |l6| . Thus short-range dis- 
order, causing intervalley scattering, needs to be present 
for AL to occur In GNRs, scattering by imper- 

fect edges may be one of the reasons behind this |8l4ll| . 
Additionally, Raman spectroscopy measurements of bulk 
graphene have revealed the presence of resonant scatter- 
ing [l^, which is another potential source of AL. The 
scattering may be due to hydrogen atoms [1^ or hydro- 
carbons 
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In this Letter, we present numerical simulations of both 
Ohmic and localized systems, showing that point-like 
scatterers can effectively be described through a formal- 
ism based on defining scattering cross sections a{E) for 



the defects. The elastic mean free path lei{E) is related 
to the scattering cross sections Ci {E) of the different de- 
fect types and the corresponding defect densities Ui via 
lei{E) = ^1 ^iniai{E) (note that in two dimensions, 
ai{E) is given in units of length). The scattering cross 
section approach provides a powerful means to estimate 
the conductance of a realistically sized GNR or carbon 
nanotube (CNT) with a finite number of point-like de- 
fects. We limit our discussion to short-range scatterers. 

For a specific defect type, the scattering cross section 
may directly be obtained from the transmission function 
T{E) of a conductor with one or several defects of the 
same type [Hjll^l- The conductance G{E) is given by the 
Landauer formula G{E) = {2e'^ / h)T {E) . In the Ohmic 
regime, the expression for <7{E) in a system with N de- 
fects is 



(j{E) = W 



TpjE) - {T{E)) 
N{T{E)) 



L « am. (1) 



where W is the width of the system (or the circumfer- 
ence of a CNT), To{E) is the transmission function of 
the corresponding pristine conductor and {T{E)) is the 
ensemble average over different defect positions and ori- 
entations. Eq. ([1]) is valid when the length L of the con- 
ductor is much shorter than the localization length £,{E). 

In the localized regime, the distribution of transmis- 
sion values is not Gaussian, but rather log-normal. Thus, 
the typical transmission Ttyp(£') = exp(logT(£^)) is a 
meaningful scaling variable. In a single-mode conductor, 
it scales as TtypiE) = exp(-L/f (£■)) Extending 
the scaling law to describe a multi-mode conductor, like 
a GNR or CNT, and treating the modes as conductors 
connected in parallel, we arrive at the expression 



TtypiE) ^To{E)cxp{-L/aE)). 



(2) 



In systems belonging to the orthogonal Wigner-Dyson 
symmetry class (e.g. graphene with short-range disorder 
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FIG. 1. (Color online) Scattering cross section cr{E) of a 
monovacancy in a (40,0)-CNT, computed from the transmis- 
sion function of a single defect (shown in the inset together 
with the transmission function of the corresponding defect- 
free conductor To{E)) using Eq. ([T|. The arrows show the 
locations of the energy values used in the statistical analysis 
of Fig. [2] The energy is given with respect to the top of the 
valence band Ev 



in the absence of a magnetic field), £,{E) is related to 
leiiE) and cr{E) through ^ 



(TojE) + l%i{E) ^ TpjE) + 1 
2 2na{E) 



(3) 



An expression for the scattering cross section, valid in 
the localized regime, is obtained by combining Eqs. ([2]) 
and ([3]): 



a{E) 



W{To{E) + l){\og{To{E)/TiE))) 
2N 



(4) 

We benchmark the formalism by comparing a single- 
defect based scattering cross section with transport re- 
sults for systems with multiple defects. An unrelaxed 
monovacancy serves as a model defect, and to exclude 
edge effects, we have chosen a (40,0)-CNT, with a circum- 
ference of about 10 nm, as a model system. Unrelaxed 
monovacancies may also be used to model adsorbate hy- 
drogen atoms (26| . as sp"^ hybridization creates a hole in 
the 77 electron network, but also more detailed models for 
hydrogen adsorbates have been used 27, 28| . 

We compute the transmission function by applying 
standard Green's function based methods to a system 
formed by a device region containing the defects and two 
semi- infinite leads |29| . The system is described by an or- 
thogonal tight-binding (TB) model, with hopping values 
obtained from Ref. [30[. A nearest-neighbor distance of 
1.42 nm is assumed, and the values are scaled to obtain 
a nearest-neighbor hopping energy of —2.7 eV, in order 
to match our previous ab initio results fs?]. All hopping 
values predicted to be smaller than 0.05 eV are set to 
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FIG. 2. (Color online) (a) and (b): Length dependence of 
the estimated scattering cross sections for a monovacancy in 
(40,0)-CNTs at the two energies shown in Fig. [1] (c) and (d): 
Sample skewnesses of the distributions of T and log(r). The 
defect density is 9.6 x 10~^ nm~^, and the vertical bars indi- 
cate the localization lengths based on the scattering cross sec- 
tion shown in Fig.[T] (e) Distribution of T(E\) at L ~ 300 nm, 
marked by the arrow in (c). (f) Distribution of log(r(i52)) at 
L « 1000 nm, marked by the arrow in (d). 



zero, which in a pristine system means that hoppings to 
farther than third nearest neighbors are excluded. The 
scattering cross section of a monovacancy in a (40,0)- 
CNT, calculated using Eq. (|T|), is plotted in Fig. (TJ to- 
gether with Tq{E) and T{E). As predicted by Fermi's 
golden rule, the van Hove singularities in the density of 
states (DOS) give rise to a greatly enhanced scattering 
rate near the band edges. 

Next, we test how the single-defect scattering cross 
section compares against results for larger systems. Figs. 
2 (a) and (b) show the estimated scattering cross sec- 
tions given by Eqs. dH) and @ as the length of the 
system is increased, keeping a constant defect density of 
n = 9.6 X 10""^ nm~^. Each point has been obtained 
from an ensemble of 5000 different realizations of de- 
fect locations. The prediction of Eq. ([T]) grows exponen- 
tially as the length of the ribbon exceeds the expected 
localization length, whereas the prediction by Eq. (|4]) 
converges towards the estimate predicted by the trans- 
mission function of a single defect. Figs. [2] (c) and (d) 
provide more insight into the behavior of the transmis- 
sion values by showing the sample skewnesses (32| of T 
and log(T). At Ei, the distribution of the transmission 
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FIG. 3. (Color online) Predicted (solid lines) and calcu- 
lated (markers) localization lengths ^{E) for (40,0)-CNTs 
with monovacancies at three defect densities. 




FIG. 4. (Color online) (a) Comparison between DFT and TB 
based conductances for a 16- AGNR with a single Stone- Wales 
defect, (b) Scattering cross sections for Stone- Wales defects, 
555777 defects and monovacancies. 



values of very short conductors containing less than ten 
defects is skewed, but even in systems with very few de- 
fects the mean value of the distribution is very close to 
the prediction based on a single defect. As shown by 
Fig. [21(e), once the length of the conductor and thus the 
number of defects increase, the distribution of T{E) be- 
comes Gaussian-like, but as £,{E) is exceeded, the skew- 
ness of the distribution starts to grow rapidly. On the 
other hand, the skewness of the distribution of log(T) 
slowly approaches a value close to zero as the the local- 
ized domain is entered. However, even at roughly ten 
times the estimated localization length, the distribution 
is still slightly skewed, and as Fig. [5] (f) shows, a rela- 
tively large fraction of the transmission values are of the 
order of 10""'^, although Ttyp is of the order of 10~^. The 
variance of the distribution shown in Fig. [21(f) equals 1.9 
times the mean value of — log(T), which is close to the 
value of two predicted by random matrix theory {25| . 

As a demonstration of the accuracy of our approach, 
we compare predicted and calculated localization lengths 
for wide ranges of energies and defect densities in Fig. [31 



The predictions are based on Eq. ([3]) and the single- 
defect scattering cross section shown in Fig. [l] whereas 
the calculated values have been obtained from systems 
ten times as long as the predicted localization lengths, 
using Eq. Each value corresponds to an ensemble 
of 200 defect realizations. At low defect density, ^{E) 
exhibits clear peaks, in accordance with Fig. [1] but at 
higher densities the peaks corresponding to the DOS of 
a pristine CNT are smoothed out. However, even at the 
highest defect density shown, the scaling approach pre- 
dicts the magnitude of £,{E) correctly. Thus we are able 
to predict properties of mesoscopic systems, only using 
the information from a single-defect calculation. 

Based on Kubo-Greenwood (KG) simulations, it has 
been proposed that topological defects, like Stone- Wales 
defects and 555777 defects, exhibit fingerprint-like scat- 
tering properties [s^. A Stone- Wales defect is a 
metastable bond rotation, consisting of two pentagons 
and two heptagons embedded in the graphene lattice, 
whereas a 555777 defect is a relaxed form of a di va- 
cancy, consisting of three pentagons and three heptagons 
3J]. Such defects are expected to be found in irra- 



diated graphene, where especially 555777 defects may 
occur in relatively high concentrations [s^. Fig. [H (a) 
compares transport results for a 16-atom wide armchair- 
edged GNR (AGNR), containing a single Stone- Wales 
defect, obtained both from a modified TB model and 
the density functional theory (DFT). The TB hopping 
parameters are obtained from the relaxed bond lengths 
around the defects, acquired from the SIESTA implemen- 
tation 3^ 37 1 of DFT and the DFT transport results are 
from the transiesta code 



In Figure [H (b), we show bulk scattering cross sec- 
tions for three defects, computed by applying /c-space 
sampled periodic boundary conditions in the transverse 
direction [40j . Although retaining the same general shape 
as the one shown in Fig. [1] the scattering cross section 
of a monovacancy is now considerably smoother. At low 
defect densities of up to roughly 10~^ defects per atomic 
site, the scattering cross section for monovacancies fits 
both qualitatively and quantitatively the corresponding 
KG based mean free paths [ilj. Also the scattering 
cross sections for Stone- Wales and 555777 defects agree 
fairly well with recent KG results although those 
are based on a somewhat differently parameterized TB 
model. When comparing against KG results, one has 
to remember that due to the phenomenon of minimum 
conductivity [l6, 42|, 2D graphene will not enter the lo- 
calized regime. At very high densities, the defects lose 
their point-like nature, and the scattering cross section 
formalism breaks down. Additionally, a high density of 
defects will contribute to the transmission through an 
impurity band [27|. 

We next test the predictive power of the scattering 
cross section approach. If a system contains several dif- 
ferent types of defects, one can estimate the average and 
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FIG. 5. (Color online) Predicted and calculated conductances 
of a realistically sized AGNR with 100 each of Stone- Wales 
defects, 555777 defects and monovacancies, as well as the con- 
ductance of a pristine conductor Tq{E). In the Ohmic (local- 
ized) regimes, rprcd(i5) refers to the estimate given by Eq. ^ 
(Eq. ([6])) and Tcaic(£) to the mean (typical) transmission. 



typical transmissions from 
To{E) 



{T{E)) 
and 



1 + n,a^{Ey 
To{E) 



exp 



2L 



To{E) + 



iT^^iME) 



L « ^{E) (5) 



L:>m- (6) 



We define the localized domain as the region where 
Ttyp{E), given by Eq. ([6]), is smaller than {T{E)), given 
by Eq. ([5]). Fig. [S] shows calculated mean and typical 
transmission functions for a 1.3 /xm long and 30 nm wide 
AGNR with 100 each of monovacancies. Stone- Wales de- 
fects and 555777 defects, together with corresponding 
predictions obtained using Eqs. ([5]) and The calcu- 
lated transmissions are based on an ensemble of 12 defect 
realizations, and the predictions on the bulk scattering 
cross sections shown in Fig. |4] (b). As the results indi- 
cate, the mean or typical transmission in a realistically 
wide GNR or CNT can be estimated by only calculat- 
ing the transmission functions for single defects in an 
edgeless system. The estimates are slightly lower than 
the actual mean and typical transmissions, which corre- 
spond to about 15% lower values for the scattering cross 
sections than predicted from systems with single defects. 
The discrepancy may be caused by the anisotropy of the 
Stone- Wales and 555777 defects, as the single-defect scat- 
tering cross section corresponds to scattering only in the 
armchair direction. 

In summary, we have presented numerical simulations 
showing that effective scattering cross sections for de- 
fects, combined with scaling theory, can be used to esti- 
mate the transport properties of graphene-based devices 



of sizes ranging from nano- to micrometers. In particular, 
their conductances can be predicted both in the Ohmic 
and strongly localized regimes. As the scattering cross 
section can be computed from small scale systems, possi- 
bilities to model systems beyond the reach of present-day 
methods are opened. 
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Academy of Finland via the FiDiPro and CoE programs. 
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